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Effects of crack tip geometry on dislocation emission and cleavage: A possible path to 

enhanced ductility 

J. Schi0tzB, L. M. Canel, and A. E. Carlsson 
Department of Physics, Washington University, St. Louis, MO 63130-^899, USA 

(June 27, 1996) 

We present a systematic study of the effect of crack blunting on subsequent crack propagation 
and dislocation emission. We show that the stress intensity factor required to propagate the crack 
is increased as the crack is blunted by up to thirteen atomic layers, but only by a relatively modest 
amount for a crack with a sharp 60° corner. The effect of the blunting is far less than would 
be expected from a smoothly blunted crack; the sharp corners preserve the stress concentration, 
reducing the effect of the blunting. However, for some material parameters blunting changes 
the preferred deformation mode from brittle cleavage to dislocation emission. In such materials, 
the absorption of preexisting dislocations by the crack tip can cause the crack tip to be locally 
arrested, causing a significant increase in the microscopic toughness of the crack tip. Continuum 
plasticity models have shown that even a moderate increase in the microscopic toughness can 
lead to an increase in the macroscopic fracture toughness of the material by several orders of 
magnitude. We thus propose an atomic-scale mechanism at the crack tip, that ultimately may 
lead to a high fracture toughness in some materials where a sharp crack would seem to be able 
to propagate in a brittle manner. 

When the crack is loaded in mode II, the load required to emit a dislocation is affected to a 
much higher degree by the blunting, in agreement with the estimates from continuum elasticity. 
In mode II the emission process is aided by a reduction of the free surface area during the 
emission process. This leads to emission at crack loadings which are lower than predicted from 
the continuum analysis of Rice. 
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I. INTRODUCTION 



The ability to calculate, or even to understand quan- 
titatively, the fracture toughness of materials has long 
been one of the major goals of materials science. Unfor- 
tunately, that goal has yet to be realized. The phenom- 
ena that cause some materials to be brittle and some to 
be ductile involve mechanisms at vastly different length 
scales, from the atomic-scale processes at crack tips and 
dislocation cores, to the formation of cell structures and 
texture at much larger length scales. 

Recently, the atomic-scale behavior of the crack tip has 
been modeled with some success. Generally, the behavior 
of cracks falls onto two categories. In "intrinsically brit- 
tle" materials, an atomically sharp crack will propagate, 
possibly leading to cleavage of the specimen. In "intrin- 
sically ductile" materials, the crack will emit one or more 
dislocations instead of propagating, and brittle fracture 
becomes impossible. Significant progress has been made 
in understanding the intrinsic behavior of the pCrack tip 
in terms of well defined atomic-scale energiesaiZl. 

The intrinsic behavior of the crack tip does not by itself 
determine the behavior of the material. Most materials 
contain a density of dislocations and dislocation sources 
that is high enough to affect fracture strongly. It is per- 
fectly possible that an "intrinsically brittle" material may 
behave in a ductile manner, if for example a loaded crack 
causes dislocations to be emitted from nearby sources, 



which subsequently shield the crack. Modeling of these 
processes is impractical at the atomic scale, but requires 
the use of continuum or quasi-continuum methodsatil. 

Even if the dislocation activity in the material sur- 
rounding a crack is in itself insufficient to change the 
behavior, the dislocations may directly interact with the 
crack tip, and thus change its behavior. If the crack 
tip intersects a dislocation a localized step or jog in the 
crack tip is formed. The jog may act as a local niide.- 
ation site for dislocations, favoring ductile behaviorE3E3. 
But even without introducing inhomogeneities the de- 
tailed atomic configuration may significantly change the 
behavior of the crack: In a recent paperlU we demon- 
strated that the emission of dislocations becomes favored 
if the crack tip is blunt at the atomic level. In this pa- 
per this effect is investigated further. This supplements 
earlier work showing that the effect of blunting on crack 
propagation is smaller than what one would expect imm 
relatively simple continuum elasticity considerationsQEj. 
We present arguments for this reduced effect based on a 
more detailed continuum model that includes the sharp 
corners of the crack. In addition, simple models of the 
inter-atomic interactions permit us to analyze the results 
in a more detailed way, and to understand the effects of 
blunting as a combination of linear elastic and nonlinear 
phenomena. 

The effect of crack blunting on crack behavior may be 
of significant practical interest in some materials. Blunt- 
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ing may occur due to emission from the tip or by collision 
of the tip with pre-existing dislocations. If this blunting 
is enough to arrest the crack the macroscopic fracture 
toughness of the material may change significantly. The 
absorption of pre-existing dislocations may be enhanced 
by attractive elastic interaction bet we ft n .the crack tip and 
the dislocations. Recent simulationsllSO have shown that 
the stress field of a crack moving through a material will 
cause dislocations located within the immediate vicinity 
of the crack path to be attracted to the crack. Since dis- 
locations lying within several hundred lattice constants 
of the crack path may collide with the tip, the likelihood 
of blunting due to intersection of the crack tip with pre- 
existing dislocations is high. 

We have also studied the behavior of the blunt crack 
under mode II loading. Here we find a much larger effect 
of the blunting on the load required to emit a dislocation, 
in agreement with arguments from continuum elasticity. 
We also find that the surface energy plays a prominent 
role in the emission process. Emission under mode II 
conditions is observed to be accompanied by a reduction 
in the total surface area. When the surface energy is 
high, emission occurs at loadings which are much lower 
than the corresponding values predicted by Riceu. 

In the next section we present the method used for 
our simulations. In section III we present the results for 
a blunt crack loaded in mode I. We find that blunting 
the crack tip results in a modest increase in the load re- 
quired to propagate it, and that often the crack response 
changes from brittle cleavage to dislocation emission. We 
discuss how this may influence the fracture toughness of 
a metal. Section |V] is a brief overview of the results in 
mode II loading, where we find that the surface energy 
plays a prominent role in the emission process. Finally, 
we present analytical calculations on stress fields around 
a blunt crack in the Appendix. 



II. SIMULATION METHODS. 

The slow decay of the elastic fields around a crack 
(a ~ r^ 1 / 2 ) makes atomic-scale simulations of cracks 
very challenging, since large system sizes are needed to 
minimize the effect of the boundary conditions. Many 
methods have been used to get around these problems, 
such as embeddiag the atomistic simulation in a con- 
tinuum model" 8 H 19 1 or using extremely large systems on 
state-of-the-art parallel computersEJ. 

Our solution is to model the elastic response of the sur- 
rounding medium by a Green's function. The simulation 
cell is divided into two zones (see figure |l|) . The atoms lo- 
cated near the crack tip interact with each other through 
a non-linear force law (the non-linear zone), whereas the 
atoms far from the crack tip interact only through linear 
forces. The displacements of the atoms in this linear zone 
can then be fully described by a lattice Green's function 
Gij(r,r'), which describes the displacement of the atom 
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FIG. 1. Division of the system into a large linear zone and a 
small nonlinear zone. Atoms in the linear zone are treated in- 
directly through the use of a Green's function. The non-linear 
zone is placed around one end of the crack. The crack is 
loaded in the middle by a force dipole, and its length L is 
kept constant at 65 lattice spacings as the blunting is varied. 



at r to a force acting on the atom at r', taking the pres- 
ence of the crack into account. The interactions of atoms 
in the non-linear zone with atoms in the linear zone are 
also described by the Green's function, and are there- 
fore linear. The Green's function can be calculated in a 
computationally efficient wayEil. The total energy of the 
system can then be described as the sum of the energy 
in the elastic far field (calculated from ifae Green's func- 
tion) and in the non-linear interactionsca. This typically 
reduces the number of degrees of freedom in the problem 
from approximately 10 7 to 10 2 -10 3 . 

We use this method to study the deformation modes of 
blunt cracks with up to thirteen atomic layers of blunting, 
for a range of force laws. Figure ||a shows a typical initial 
configuration. We have a region in front of the crack 
where the atoms are allowed to move freely. Two "spurs" , 
along which dislocations generated at the crack tip can 
move away from the crack, are included in the non-linear 
zone. Since bonds cannot be broken and reformed in 
the linear zones, dislocations will be unable to leave the 
non-linear zone. The two-dimensional hexagonal lattice 
studied contains three slip planes at 60° angles. The 
crack is oriented parallel to one of the slip planes. 

The inter-atomic interactions in the npii-linear zone 
are described by the UBER pair potentiate: 



F(r) — —k(r — ro) exp 



r - r 
I 



(1) 



where r is the inter-atomic separation, ro is the sepa- 
ration in equilibrium, and I is a range parameter. The 
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interactions are cut off so that only nearest-neighbor in- 
teractions are included, and forces are shifted slightly to 
avoid steps in the force at the cutoff. Further, a small 
scaling C of the force law is used to preserve the elastic 
constants, thus enabling us to use the same Green's func- 
tion for all the force laws. The force law thus becomes 



F{r) = C 



—k(r — rp) exp 



r - r 



F 



ro 



r C ut ~ r 



(2) 



where Fq = —k(r cut — ro) exp ((r cut — »"o)/0 assures that 
the force is zero at the cutoff distance r cut (in this work 
1.7tq). Since we do not study specific materials, no 
attempt is made to use realistic many-body potentials. 
In spite of the relatively simple potentials used, this 
methodology is known to give results that agree well with 
the observed behavior of real materials, and with contin- 
uum modelsdO. We also see good agreement between 
our calculaticoisi-and calculations using the E mbedded 



Atom MethodJa, as discussed in section [II B 



For all the calculations presented in this paper, the to- 
tal system size has been 1024 x 1024 atoms. In order to 
establish that the results are not influenced by the bound- 
ary conditions, we have repeated a few of the simulations 
with 16 times as many atoms (4096 x 4096). The only 
difference is a consistent 1-1.5% reduction of the loads 
required to cleave the crack or to emit a dislocation. All 
trends remain the same. The size-dependence is caused 
by the stress fields of the cracks in neighboring supercells 
causing the true stress intensity factor to be slightly dif- 
ferent from the calculated values. Since this interaction 
remains the same it does not influence our results, where 
we study the changes in the critical loads. We have also 
verified that increasing the size of the non-linear zone 
does not affect the results. 



III. THE BLUNT CRACK UNDER MODE I 
LOADING. 

A. Simulations. 

When we load the cracks until an instability occurs, 
one of two different responses are observed: cleavage and 
dislocation emission (figure |J) . The dislocation emission 
always occurs in the downwards direction, as shown in 
the figure, except for the case of a single layer of blunt- 
ing, where the crack geometry is symmetric. Further- 
more, the emission preserves the shape of the blunted 
crack. This means that further dislocations can be ex- 
pected to be emitted in the same direction, provided that 
the dislocation can move so far away that its effect on the 
local stress field at the crack tip is small. The crack shape 
studied is thus preserved during the blunting process. 

Figure ^| shows the force required for cleavage or dis- 
location emission as a function of the blunting, for seven 
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FIG. 2. The two instability modes of the crack, (a) shows 
the configuration before plastic deformation. The crack, three 
atomic layers thick, extends out of the left of the picture. 
Only atoms in the non-linear zone are shown, see text, (b) 
shows the crack propagating by cleavage, (c) shows the crack 
emitting a dislocation. The dislocation has traveled along the 
black line. Since it cannot leave the non-linear zone, it has 
been pinned at the bottom (indicated by the arrow) . To make 
the path of the dislocation visible, the atoms were given two 
different shades. Prior to emission the atoms were lined up in 
rows of the same shading. Where the dislocation has traveled 
lines of atoms of different shade meet. 



different decay lengths I of the force law. The results for 
the sharp crack (zero blunting) are reported on the y- 
axis. We find that increasing the decay length I increases 
the load required to propagate the crack (see also Zhou 
et alB). This is because increasing I raises the surface 
energy and therefore the loading to cleave. The critical 
values observed are 4-7% above the values expected from 
the Griffith criterion, this discrepancy is attributed to the 
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FIG. 3. The external load required for cleavage or disloca- 
tion emission as a function of the blunting, for seven different 
force laws. Open symbols indicate dislocation emission, closed 
symbols indicate cracking. It is seen that for a wide range of 
force law parameters the sharp crack propagate by cleavage, 
whereas the blunt crack emits dislocations. For Z > 0.5 the 
crystal structure ahead of the crack changed from triangular 
to square under influence of the stress, making further simu- 
lations meaningless. 



discreteness of the lattice (lattice trappingo) . Two main 
effects are seen with increasing blunting. For all force 
laws, blunting increases the critical loadings to cleave or 
to emit, but only by a modest amount. For 13 layers 
of blunting the increase is only approximately 20%. This 
should be compared to estimates derived from linear elas- 
ticity of elliptic cracks giving an effect of around 40% for 
a single layer of blunting, and 250% for 13 layers!! Sec- 
ondly, for a range of force law parameters / a normally 
brittle crack is observed to emit dislocations after being 
blunted by one or two layers. Here, although the force 
required to cause the blunt crack to emit is still larger 
than the load that cleaves the sharp crack, the force re- 
quired to propagate the crack has increased even more, 
shifting the balance between emission and cleavage. 

For each value of the blunting, there is a transition 
from cleavage to dislocation emission with increasing L 
caused by an increase in the surface energy. Figure |] 
shows the critical value of the force law parameter l c for 
which the crack response changes from cleavage to emis- 
sion. The first layer of blunting results in a significant 
reduction of l c , causing a wide range of force laws to fa- 
vor emission. Further increase in the blunting causes l c to 
increase marginally. Thus the first layer or two of blunt- 
ing tend to enhance ductility, and the subsequent effects 
are much smaller. We have not been able to explain this 
effect on the basis of analytic theory. 
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FIG. 4. The force law range parameter l c where cleavage 
and dislocation emission are in equilibrium, as a function of 
the blunting. The first layer of blunting has a large effect, 
making a range of force laws (l c > 0.27) emitting, further 
blunting makes some of these cleave again. 
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FIG. 5. The exponent of the stress singularity as a func- 
tion of the opening angle for three different loadings. See the 
Appendix. 



B. Discussion of the mode I results 

In order to understand these results, we have analyzed 
the nature of the elastic fields around a 60° wedge crack 
(see the Appendix) . Simplistically, one would expect that 
the increased opening angle should reduce the stress sin- 
gularity at the crack tip, making cleavage more difficult. 
However, we find that the strength of the stress singu- 
larity at the end of a blunt crack with a 60° angle is not 
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FIG. 6. The components of the elastic strain five atoms 
ahead of the crack as a function of blunting. The loading of 
the crack is kept constant. 
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FIG. 7. The effect of reducing the non-linear zone at the 
crack tip. When the atoms shown in black are removed from 
the non-linear zone (i.e. they are restricted to interact through 
linear forces), the effect of the blunting is reduced. The range 
parameter in the force law is I — 0.25. 



changed significantly compared to a sharp crack. In par- 
ticular, the exponent of the singularity has only changed 
from -0.5 to -0.488, see figure S. Therefore, purely from 
considerations of the crack tip field, we would not expect 
any significant change in the load required to propagate 
the crack, and believe that this explains the smallness of 
the blunting effect. This expectation is confirmed by fig- 
ure ||, which show the elastic strain field five atoms ahead 
of the crack for a completely linear (but still discrete) sys- 
tem. The effect of the blunting on the tensile strain (e yy ), 
which should be the most important for cleavage, is mini- 
mal, although some reduction is seen when the first layer 
of blunting is introduced. When the crack tip is blunted 
only by one layer, the crack is symmetric and it does not 
resemble the shapes studied in the Appendix; the behav- 
ior near a 60° corner is thus not relevant, nor is contin- 
uum theory expected to be a good description when only 
a single atomic layer has been removed. As the blunting 
is increased, the actual crack shape approaches that of 
the ideal theoretical shape studied in the Appendix, and 
the stress field then approaches the value corresponding 
to the wedge crack. Since that value is virtually indistigu- 
ishable from the stress in front of a sharp crack, the stress 
is seen to increase again when the blunting increases. 

The decrease in strain when the first layer of blunt- 
ing is introduced is consistent with the observed increase 
in the force required to propagate the crack. However, 
the elastic behavior completely fails to explain the ad- 
ditional increase in the force to cleave as the crack is 
further blunted, since it would—indicate that this force 
should decrease slightly againnJ. Clearly other effects 



are important. 

One such effect is the non-linearities in the force law. 
We have found that these nonlinearities explain the in- 
crease in the cleavage stress intensity seen in figure [| To 
evaluate their effects, we have partly "turned them off" 
by reducing the size of the non-linear zone. The results 
are illustrated in figure for I = 0.25. Here the non- 
linear bonds (cf. inset) include only the bonds between 
surface atoms. Reducing the size of the non-linear zone 
does not change the critical loading when the blunting is 
only one atomic layer. For larger values of the blunting, 
reducing the non-linear zone reduces the critical load to 
cleave. The effect of the blunting the crack beyond one 
layer thus becomes smaller when the non-linear zone is re- 
duced. This is caused by a stretching of the bonds of the 
atoms near the end of the crack. The non-linear force-law 
permits bonds to be stretched further than would other- 
wise be the case, thus accommodating more of the strain. 
This will tend to unload the crack tip bond, inhibiting 
cleavage. 

This leads us to conclude that the inhibition of crack- 
ing caused by the first layer of blunting is mainly due to 
a change in the elastic fields around the crack, but that 
further inhibition is caused by non-linear bonds along the 
end of the crack unloading the crack tip. 

Since the effects of the nonlinearities can be expected 
to depend somewhat on the actual potential used in the 
simulations, it is not entirely clear that these effects will 
be present in three-dimensional simulations using more 
complicated potentials. The effect of the changes in the 
elastic field are going to be preserved, so we can expect 
that the effect of the first layer of blunting will be pre- 
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served. However, it is unclear how effects like surface 
stress and surface relaxations, typically seen both in ex- 
periments and in many-body potentials, will change the 
behavior as the blunting is increased. Here it is reassuring 
to note that recent simulations of cracks in NiAlHa and 
in A13 give results very similar to the ones observed here. 
Both sets of simulations use the Embedded Atom many- 
body potentials. The crack geometries differ somewhat 
from this work, but the increases in the critical loads are 
similar. For one crack orientation the shift from cleavage 
to dislocation emission is observed in NiAl, similarly to 
the results presented here. 

It is now clear that the presence of sharp corners in the 
blunt cracks is highly significant for its behavior, since 
stress concentrations appear at these corners that may 
be almost as strong as the stress concentration at the tip 
of a sharp crack. Macroscopically blunted cracks may 
appear to have a smooth shape, but when the blunting 
is only a few tens of lattice constants, or even smaller, 
the crystal structure makes it almost inevitable for sharp 
corners to be present in the crack tip configuration. Fur- 
thermore the dislocation emission processes creating the 
blunting will almost unavoidably lead to angular crack 
tip shapes similar to the one studied here, and will not 
create smoothly blunted cracks. The exact angles of the 
corners will, of course, depend on the detailed crystal- 
lography. From figure || it can be inferred that for 90° 
and 120° corners, the effect of the blunting on the stress 
fields is still small, although the effect is expected to be 
significant at least in the 120° case. 

The change in deformation mode, from a brittle cleav- 
age to a ductile emission as the crack is blunted on the 
atomic scale, may have a significant impart, cm .the macro- 
scopic behavior of the material. MesarovicESEZI has shown 
how a crack tip moving through a material will attract 
dislocations located in a relatively broad strip ahead of 
the crack tip, causing them to move towards the tip and 
ultimately collide with it. Such collisions will blunt the 
crack, and in some materials be sufficient to arrest it. 

The two-dimensional nature of the simulations pre- 
sented here invoke the picture of infinitely long straight 
dislocations interacting with a straight crack tip. In real- 
ity neither the crack nor the dislocations have such ideal- 
ized shapes. The dislocations interacting with the crack 
are likely to be short segments or loops, and only a short 
segment of the crack is thus likely to be blunted. The 
result is a crack front where short segments are arrested 
while the rest of the crack moves on, until either the lo- 
cal stress intensity factor becomes sufficient to reinitiate 
crack propagation, or until the sharp crack moves around 
the obstacle. This local "crack trapping effect" gives rise 
to a gradually increasing ceack tip toughness as the dislo- 
cation density is increasedlHl, instead of the complete in- 
hibition of sharp crack propagation that two-dimensional 
simulations would suggest. 

Recently, Beltz et alc3 have proposed a continuum 
plasticity model for cleavage. They show that the "shield- 
ing ratio" , i.e. the ratio between the macroscopic tough- 
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FIG. 8. Left: Geometry of the non-linear zone under mode 
II loading. Right: Due to the asymmetry of the crack, it is 
necessary to keep track of the sign of the applied mode II 
load. The loading indicated in the upper figure will be called 
positive, the loading in the lower figure negative. The figures 
also indicate how dislocations are emitted in the two cases, 
see text. 



ness and the microscopic crack tip toughness, depend 
sensitively on the ratio between the microscopic crack tip 
toughness and the yield stress. This means that a mod- 
est increase in the microscopic toughness of the crack tip 
may give an enormous increase in the macroscopic tough- 
ness of the material. MesarovicEZI uses this model to show 
that if collisions with pre-existing dislocations locally ar- 
rest the crack, they can increase the crack tip toughness 
by a factor 5-7 causing an increase of the macroscopic 
toughness by 2-3 orders of magnitude. In this way he 
can explain the vast difference in fracture toughness of 
iron-silicon and tungsten at 77K. 

This mechanism, where collisions with preexisting dis- 
locations blunt the crack and thereby arrest segments of 
the crack, will only lead to a high macroscopic toughness 
if the blunting arrests the crack completely, or signifi- 
cantly increases the load required to propagate the crack 
further. We find that the complete arrest is a possibility, 
since cracks in some materials (those with high values of 
I in this study) will start emitting dislocations in stead of 
propagating as soon as the crack is blunted. Other mate- 
rials, however, will not change behavior. In these materi- 
als the effect of blunting is negligible, and the mechanism 
discussed above will not lead to a high toughness. 



IV. THE BLUNT CRACK UNDER MODE II 
LOADING. 

We have also investigated the behavior of blunt cracks 
under mode II loading. Figure ^ shows the non-linear 
zone used for these simulations. In this case the spurs al- 
lowing dislocation emission on inclined planes have been 
removed. Since the crack is not symmetric, the sign of 
the load is important. The sign convention used is given 
in figure ||. When the crack is loaded, it is seen that 
dislocations are not necessarily emitted in the sharpest 
corner. Under negative load the dislocation is emitted 
in the lower, more blunt corner, where the stress concen- 
tration is smaller (see figure . This is an effect of the 
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FIG. 9. (a) Dislocations emitted from the upper corner un- 
der positive load, (b) Dislocations emitted from the lower cor- 
ner under negative load, (c) Collapse of the corner of the crack 
under positive mode II loading for sufficiently high blunting, 
(d) Emission of a dislocation at a slip plane above the ex- 
pected plane, accompanied by a reconstruction of the corner. 
Note: In the two upper figures the initial blunting was three 
layers, in the lower figures it was four layers. In all the fig- 
ures the atoms have been given two different colors, forming 
straight lines before the dislocation was emitted. The path 
followed by the dislocation can thus easily be seen. 



surface energy, in both cases the dislocation is emitted 
in such a way that the surface area is reduced by the 
emission. The surface energy thus assists the dislocation 
nucleation, and this effect dominates that of the stress 
singularity. This effect is similar to the "ledge effect" 
under mode I loading, where the creation of extra sur- 
face at the crack tip may dominate the energetics of the 
emissioncl. 

Figure [h] shows the load required for mode II emis- 
sion, under positive loading. No measurements were at- 
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FIG. 10. The external load required to emit a dislocation 
under positive mode II loading, for a variety of force laws. 
Symbols connected with a solid line are data points obtained 
with the full nonlinear zone. Symbols connected by a dashed 
line are data points calculated with a reduced nonlinear zone. 
Open symbols indicate a slightly different emission mode, see 
text. For most force laws, the crack collapsed if the blunt- 
ing was only one layer, and no data could be obtained. The 
symbols on the y-axis are the Rice prediction for the sharp 
crack. 



tempted for the sharp crack, since it would require a si- 
multaneous mode I loading to keep the crack open. In 
the graph, filled symbols connected by a solid line indi- 
cated data obtained using the full non-linear zone shown 
in figure |^. For four or more layers of blunting the-corner 
collapsed when a load was applied, see figure f*£j. The 
collapse could be prevented by using a smaller nonlinear 
zone, where the upper and lower two layers in figure |8| 
are removed from the nonlinear zone. The data obtained 
in this way are shown connected by a dashed line. The 
reduction of the non-linear zone is seen to cause a tiny 
change in the force required to emit a dislocation. For 
the force law I = 0.25 the corner collapsed for all levels of 
blunting, only data obtained with the reduced nonlinear 
zone is given. For I — 0.30, dislocation emission again 
occurred before collapse of the corner for sufficiently low 
levels of blunting. For all force laws with I > 0.25 the 
emitted dislocation immediately climbed one slip plane 
up, leaving an atom in the corner (figure ^d). These cases 
are shown with open symbols on the graph. They clearly 
have a different dependence on the blunting, and since 
the dislocation moves along the very edge of the reduced 
non-linear zone, such reduction has a significant effect on 
the force required to emit. 

The results of loading the crack in the negative direc- 
tion is shown in figure [0]. Since the material is very 
brittle for I = 0.15, the crack preferred to branch and no 
emission data were obtained for that force law. 
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Blunting (layers) 

FIG. 11. The external load required to emit a dislocation 
under negative mode II loading, for a variety of force laws. 
For most force laws, the crack collapsed if the blunting was 
only one layer, and no data could be obtained. No data was 
obtained for the / = 0.15 force law, since the crack branched 
instead of emitting a dislocation. The symbols on the y-axis 
are the Rice prediction for the sharp crack. 



For both directions of the load, it is seen that the effect 
of the blunting; is far more pronounced than in the mode 
I case (figure 0). This is in agreement with the expecta- 
tions from continuum linear elasticity, as shown in figure 
[5]. The stress singularity near the 60° corner is reduced 
from r -0 - 5 to r~ - 27 , giving a significant reduction of the 
stress at the corner as the blunting is increased. For the 
force laws / = 0.15 and I — 0.20, the magnitude of the 
effect of the blunting is consistent with this change in 
the exponent, if we assume that the dislocation is emit- 
ted when the shear stress one lattice constant ahead of 
the crack tip reaches a critical value. For the other force 
laws the unusual emission process mentioned above gives 
a different dependence of the blunting. 

For negative loading, emission is seen from the 120° 
corner. The dependence on blunting is even stronger 
than for positive loading. In this case continuous lin- 
ear elasticity predicts no stress singularity in the corner, 
and indeed no stress at the corner, provided the load is 
pure mode II. However, the symmetry axes of the corners 
are rotated with respect to the load, so in neither case is 
the load pure mode II as assumed in the Appendix, but 
is a mixture of mode I and mode II. 

For both signs of loading, the effect of the sur- 
face energy-, is clearly seen by comparing with the Rice 
predictionsu for the sharp cracks—as. modified to reflect 
the plane stress situation studiedoEj (the original crite- 
rion was developed for plane strain) 

K H e = V27« S M1 + v) (3) 



where a us is the so-called "unstable stacking fault 
energy"H In figures [h] and [n] the Rice predictions are 
plotted along the y-axes. These predictions have been 
shown to be in igood agreement with atomic simulations 
of sharp cracksE-3, and predict only a small variation of 
the critical load for emission when the force law is varied, 
since j us does not vary much for the values of I studied 
here. The far larger observed variation of the critical 
load is caused by the surface energy, and it is seen that a 
high surface energy (high value of I) enhances the emis- 
sion significantly. The effect of the surface energy on 
the emission process will depend on the ratio between 
the surface energy and the unstable stacking fault energy 
(is /jus), since j us controls the energy required to create 
the dislocation. For most fee metals it is approximately 
5-10, for bec metals l-7u; in our simulation it varies from 
1.4 (/ = 0.15) to 4.3 (I = 0.50). We would thus expect 
that the effect of the surface energy will be at least as 
significant in real materials as in our experiments. 

As is seen in figure |^, the stress singularity in mode III 
is stronger that in mode II, but weaker than in mode I. 
We would therefore expect that a blunt crack loaded in 
mode III behaves in a way that is intermediate between 
the observed behaviors of mode I and mode II cracks, 
i.e. we would expect a dependence on the critical load 
on the blunting that is stronger than for mode I, but 
weaker than for mode II. Since the surface area is not 
changed under mode III emission, we would not expect 
the strong dependence on the surface energy that we ob- 
serve in mode II. 



V. CONCLUSION. 

We have examined the behavior of blunted cracks un- 
der mode I and mode II loading. The mode I results in- 
dicate a new mechanism for enhanced ductility in some 
materials. We find that blunting the crack increases the 
force required to propagate the crack, but only by a sur- 
prisingly modest amount. Ten layers of blunting causes 
an increase of 15-20%. The relative modest magnitude 
of the effect can be attributed to the sharp corner of the 
blunt crack, where the stress singularity under mode I 
loading is almost as strong as for a sharp crack. Despite 
the small effect on the cleavage criterion, for many force 
laws the blunting causes the crack to change behavior 
and to start emitting dislocations instead of propagat- 
ing. This has the consequence that if a crack in such a 
material absorbs a dislocation at the crack tip, a segment 
of the crack will locally be arrested, leading to an increase 
in the microscopic crack tip toughness, which further may 
lead to a much larger increase in the macroscopic tough- 
ness of the material. The blunting effects presented here 
may thus, in combination with the effect of attracting 
dislocations to the crack tip described by MesarovicE3EZl, 
cause an increase in the fracture toughness of materials 
by many orders of magnitude. 
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Under mode II loading dislocation emission always oc- 
curs in such a way as to decrease the total surface area of 
the crack. The surface energy will thus aid the emission 
process, and may dominate over the stress singularity. 
The emission is thus not necessarily from the corner with 
the larger stress concentration. When the surface energy 
is high, the emission is seen to occur at loads significantly 
below the value predicted by Rice. The effects of blunt- 
ing are greater than for mode I loading, in agreement 
with the predictions of linear elasticity. 




FIG. 12. Left: A wedge crack with opening angle a. The 
grey area is the material. Right: The geometry of a blunt 
crack. 
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APPENDIX: ANALYTIC TREATMENT OF THE 
BLUNT CRACK. 

In order to better understand the simulations, we have 
investigated the behavior of the stress field near a blunt 
crack, in the framework of linear elasticity. Unfortu- 
nately, an exact analytic solution for the stress field 
around a blunt crack loaded in mode I or II can. only be 
found for an elliptic crack in an infinite mediumEil. In this 
section, we show that for some geometries the problem 
can be solved under the assumption of antiplane strain 
(i.e. mode III loading). Wc then use simpler geometries 
to elucidate how the change to plane elasticity (mode I 
and II) will modify the results. 



1. The blunt crack in antiplane strain. 

The stress field around a wedge-shaped crack loaded in 
mode III (see fig |l2|a) is well knownE3E3. The stress field 
displays a singularity at the crack tip, but the singularity 
is weaker than for the sharp crack, and vanishes as the 
opening angle a reaches ir. 



P 



o 



2tt — a 



(Al) 



The real geometry of a blunt crack is more like the one 
shown in figure [l^ j. Here we would expect the field far 
from the tip to be unaffected by the blunting, i.e. to scale 
like a ~ r" 1 / 2 . Near the tip the field must be similar to 
the field from a wed ge crack, i.e. a ~ r p , where p is 
given in equation Al. A weaker singularity will also be 



present at the lower corner of the crack. We can then 
approximate the total stress field by a field made up of 
these parts, by matching the stresses from the different 
solutions at an "appropriate" distance from the crack tip. 
It can be expected that this matching distance will be of 
the order of d. We will verify this naive picture with an 
exact solution of the problem, using conformal mappings. 

In antiplane strain, the stresses can be written as the 
derivatives of an analytic function, u{z) = a yz + ia xz = 
2r]'(z) where z = x + iy. This representation has the 
nice property, that the boundary conditions (traction free 
surfaces) are preserved under a conformal mapping. If 
we know fj(£) is a solution to a known problem (f2^), and 
if the function z — z(£) maps the known problem into a 
new problem (Q z ) in such a way that the boundary maps 
onto the boundary, and the material region maps onto the 
material region, then T](z) = f)(£(z)) is a solution to the 
new problem (provided that z(£) is analytic everywhere 
inside the material)! 3 ^!. 

We will choose the upper half plane as our known prob- 
lem, and use the Schwarz-Christoffel (S-C) transforma- 
tion to generate the map. Taking the angles as 7r/3 and 
27r/3 a possible S-C transformation is 



z(0 = A J dte /3 (Z-l) 1/3 + B 



(A2) 



where A and B are complex constants to be determined 
later to give the right size, orientation and position of the 

1 /3 

blunt crack. Using the transformation t = ((£— , 
and decomposing the resultant rational function into par- 
tial fractions, we get 



7T ' 



k Log(i — fc) 



1 



fe£{l,a,a} _ 



t : 



2(t- ky 



e-1 



1/3 



, (A3a) 



(A3b) 



where we have set A and B to 9/7T and 0, respectively, 
to get a blunt crack of width 1 with the crack tip at the 
origin. The constants a and a are equal to exp (±|7ri), 
i.e. the fc-sum is over the roots of fc 3 = 1. 

In f2^ (the upper half plane) the following solutions 
are consistent with the boundary condition, traction free 
surfaces: 
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Vn(0=CnC, neZ\{0}, 



Cn e C. 



(A4) 



Each gives a solution to the blunt crack, rj n (z) = 
fj n (t;(z)). However, we must discard all solutions with 
n > 2 because the stress would diverge at infinity, and 
all solutions with n < — 1 because the displacement di- 
verges at the origin. This leaves us n = 1 (n = is not 
an option, since it gives a constant i] and thus no stresses 
anywhere). 

We will now expand the solution in the near and far 
field approximations. In both cases it is important to 
keep the branch cuts of the logarithms and the non- 
integral powers out of the material region, and to avoid 
moving between different branches of these complex func- 
tions. 

For the far field, i.e. for z — > oo and £ — > oo we get 



2(0 



9 

'2?r 



£(*) 



/2vr 



and the elastic field becomes 



a(z) = 



K in 1 
%/2vr \Tz 



(A5) 



(A6) 



The far field solution is thus identical to the stress field 
around a sharp crack. From this expansion we get the 
constant c = —^Kttt. 

The field near the crack tip (z — > and £ — > 0) can 
be found by expanding equation ( |A3| ) to fifth order in 
(— z) 1 / 3 (It vanishes to fourth order!). We then get 



*(0 

and get the stress field 



p -37Ti/5 f ^ 

1 27 



3/5 



,3/5 



a(z) 



- „-™/10 



Kmjz 



-2/5 



7 = 



(57r) 2 /5 



(A7a) 
(A7b) 

(A8) 



We see that we find the same behavior as near a wedge 
crack with the same opening angle. 
The exact solution is 



a(z) 



-i-KmZ\z) 
ir 



(A9) 



whe re £,'(z) must be found numerically, since equation 
( [A3] ) cannot be inverted analytically. Figure [l3| compares 
the exact solution to the asymptotic expressions; it is 
seen that the asymptotic approximations are quite good 
everywhere on the x axis. 



2. The blunt crack in plane elasticity. 

The field equations of plane elasticity (plane stress and 
plane strain) do not allow the same elegant treatment 
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FIG. 13. The stress field ahead of a blunt crack, loaded in 
mode III. 



as in antiplane strain. Although a conformal mapping 
technique is possible^, the boundary conditions are not 
preserved under the mapping transformation, complicat- 
ing the process considerably. In many cases such an ap- 
proach can still be fruitful, especially .if .combined with 
the integral method of MuskhclishviliEJS, but the sharp 
corners of the blunt crack causes unwieldy singularities 
to appear. If the conformal mapping is expanded as a 
(truncated) power series, the problem of a blunt crack is 
soluble using Muskhelishvili's integral methocEj, but the 
resultant series expression of the stress fields can gener- 
ally not be summed analytically. 

Whereas the blunt crack cannot be solved exactly, a 
solution for the wedge crack can be found, using a differ- 
ent approach. We describe jtiie-fitress state of the system 
by its Airy stress functionEaO $. The stress function 
must satisfy the biharmonic equation V 4< & = 0, and the 
stresses are then given as suitable derivatives of $. 

We are interested in solutions to the biharmonic equa- 
tion that fall off as a power law. In polar coordinates 
such solutions have the form 



$(r,0) = r x+1 f(6) 

where 

/a (6») = a sin(A + 1)6» + c 2 cos(A + 1)6* 
+c 3 sin(A - 1)9 + c 4 cos(A - 1)6* 

giving the stresses 

1 9$ 1 d 2 $ 

r dr r 2 dQ 2 
= r A - 1 (/i / (0) + (A + l)/ A (0)) I 



(A10) 



(All) 



(A12a) 
(A12b) 
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dr \r 86 



»A-1 \ f i 



(Al2c) 



We now demand traction-free surfaces, i.e. ogg = a r g = 
for 9 = ±/3 where is given by the opening angle 
(figure |l|): (3 = n-a/2. We thus demand that f{±0) = 
f(±0) = 0, giving 

c 2 cos((A + l)(3) + c 4 cos((A - 1)0) = 0, (Al3a) 



c 2 (A + 1) sin((A + 1)0) + c 4 (A - 1) sin((A - 1)0) = 0, 

(A13b) 

d sin((A + 1)) + c 3 sin((A - 1)0) = 0, (A13c) 

ci(A + 1) cos((A + 1)0) + c 3 (A - 1) cos((A - 1)0) = 0. 

(A13d) 

These two sets of two homogeneous linear equations 
only have non-trivial solutions if one of the determinants 
vanish. For the first two equations this happens for 



Asin(2/3) +sin(2A/3) = 
and we get the stress fields 

-^^(-(A + ricos^A + l)^ 



(A14) 



c 2 r 



+ (A-3)Q A cos[(A-l) 
&66 = c 2 r A_1 A(A + 1) [cos [(A + 1)0] 



(Al5a) 



a r g = c 2 r 



where 



-Q A cos[(\-l)6]\ (A15b) 
^((A + ^sin^A + l)^ 

- (A - l)Q J 4sin[(A - 1)9] 



(Al5c) 



Q 



c 4 _ cos [(A + 1)0] 
~ 2 ~ cos[(A-l)/3] 



The solutions of equation ( A14 ) are constrained by the 
requirement that the stresses vanish at infinity and that 
the displacements remain finite at the crack tip: < A < 
1. For the sharp crack (a = and (3 = tt) we get A = 1/2 
and find Ihe known stress field of the crack under mode 
I loadings. 



Equations (Al3c) and (Al3d) give another eigenequa- 
tion: 



Asin(2/3) - sin(2A/3) = 
giving the stresses 



(A17) 



a rr = cir A " 1 A(-(A + 1) sin[(A + 1)6] 

+ (A-3)Q B sin[(A-l)0]V (Al8a) 
agg = cir A - 1 A(A + l)(sin[(A + l)6»] 

— Qb sin[(A — 1)6] ^ , (A18b) 
a r g = -cir x ~ l \U\ + 1) cos[(A + 1)6] 

+ (A- l)Q s cos[(A- 1)0]), (Al8c) 



where 



Qi 



eg = sin[(A + l)/3] 

Cl 



[(A-l)/3] 



sm 



(A19) 



Again A is constrained to < A < 1, which can only be 
fulfilled for opening angles less than a < 103°. For the 
sharp crack we recover the mode II solution. 

Figure H shows the exponent (A — 1) of the stress singu- 
larity as a function of the opening angle for wedge cracks 
loaded in mode I, II and III. 

The behavior of the mode II solution at a = 103° is 
peculiar. As the angle approaches 180° the mode I solu- 
tion approaches a half-plane under uniform tensile stress. 
Such a behavior is not possible in mode II, since the 
boundary conditions prevent a uniformly sheared half- 
plane (an attempt to apply shear alone would only result 
in a rotation). It is thus clear that the stress must go 
to zero as the opening angle approaches 180°. However 
that behavior is seen already at 103°. Between 103° and 
180° a stress field that decays as a power law in r is 
not possible, nor is a stress constant in r. It is hard to 
imagine another form of the stress where it increases as 
r decreases, since any other form would require that an- 
other length scale enters, and linear elasticity is known 
to be scale- invariant. 

As a curiosity it can be mentioned that this peculiar 
behavior also has been seen in a related problem: In the 
beginning of the century the problem of a wedge lc 
at its apex by a couple of moment M was solvecOl 
A solution was found with the resulting stress field de- 
(A16) cavm g as r~ 2 , diverging for (in our notation) a = 103°. 

This strange singular behavior at a physically meaningful 
angle was later investigated by Sternberg and KoiterCJ. 
They replaced the couple by a physically realizable force 
distribution of size o, and then solved the stress field for 
r ^> a. They found that for a > 103° the solution con- 
verged to the original solution, but below 103° an extra 
(dominating) term appeared, making the concept of a 
couple loading unphysical. That extra term is identical 
to the mode II solution above! 
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